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A note on Newtonian, Lagrangian, and Hamiltonian 
dynamical systems in Riemannian manifolds. 

R. A. Sharipov 

Abstract. Newtonian, Lagrangian, and Hamiltonian dynamical systems are well 
formalized mathematically. They give rise to geometric structures describing motion 
of a point in smooth manifolds. Riemannian metric is a different geometric structure 
formalizing concepts of length and angle. The interplay of Riemannian metric and 
its metric connection with mechanical structures produces some features which are 
absent in the case of general (non-Riemannian) manifolds. The aim of present paper 
is to discuss these features and develop special language for describing Newtonian, 
Lagrangian, and Hamiltonian dynamical systems in Riemannian manifolds. 



1. Force field of Newtonian dynamical system. 

The primary and most transparent way of describing real mechanical systems 
is based on Newton laws. Newton's second law yields differential equation for the 
motion of small particle with mass m under the action of force F: 



>■ ! TO-i^ = F(r,r). (1.1) 

^vj . Here r — r{t) is a vector of three-dimensional geometric space marking position 

[^^ ' of moving particle. Formally, one can consider the equation (1.1) for r e R" and 

^^ , can take m = 1 for the sake of simplicity. Further one can replace M" by arbitrary 

'"^ ' smooth manifold M and write the equation (1.1) in local coordinates x^, ... , a:": 

^- i'' :^^''{x\... ,x",i\... ,x"), A:=l, ...,71. (1.2) 

Cj ' Once the equations (1.2) are written, we meet the problem if interpreting these 

.. ', equations. If x^, . . . , a;" are coordinates of moving point p = p{t) in M, then their 

K^ ■ first derivatives i^, ... , i" are components of velocity vector v G Tp{M). But 

k>( I second derivatives x^, ... , x" are not components of a tangent vector of Tp{M). 

; I ' Therefore we are to consider the pair q = (p, v) being a point of tangent bundle 

C^ I TM, and then write the equations (1.2) as a system of first order ODE's: 






(1.3) 



Ordinary differential equations (1.3) correspond to the following vector field in TM: 

* = v' • A + • ■ • + ^^" • ^ + *' • A + • • • + *" • ^- (1-4) 

dx^ dx"^ dv^ dv" 
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li q = ija, v) is a point of tangent bundle TM and if tt : TM ^ M is a map of 
canonical projection, then, applying associated linear map tt* : Tq{TM) — > Tp{M) 
to the above vector (1.4), we obtain the equality 

7r»$ = V. (1.5) 

Definition 1.1. Vector field $ in tangent bundle TM satisfying the condition 
(1.5) is called Newtonian vector field. 

Definition 1.2. Newtonian dynamical system in smooth manifold M is a dynami- 
cal system determined by some Newtonian vector field in TM . 

For the motion of real particle both vectors v and F are in the same space. We 
can measure their lengths and the angle between them. Passing to general case of 
n-dimensional smooth manifold M, we loose this opportunity. Indeed, vector $ is 
2n-dimensional vector tangent to TM ^ while v is n-dimensional vector tangent to 
M . This situation changes crucially if we take Riemannian manifold M . In this 
case we can consider vector F with components 

n n 

^fe^^fe^^^pfc ^«^j^ A: = l, ...,n. (1.6) 

It is tangent to M at the point p = 7r((j). But its components (1.6) are functions 
of double set of arguments x^, . . . , a;", f ^, . . . , w". In other words, F is a vector 
in Tp(M) depending on the point q = (p, v) G TM . 

Definition 1.3. Extended vector field F in M is a vector-valued function that 
maps each point q ^ G C TM to a vector of tangent space Tp{M), where p = Tr(q). 
Subset G of TM is a domain of extended vector field F. If G = TM, then F is 
called global extended vector field in M. 

Vector F with components (1.6) is called force vector. It determines force field of 
Newtonian dynamical system in Riemannian manifold. Force field F of Newtonian 
dynamical system is an extended vector field in the sense of definition 1.3. Velocity 
vector V can also be treated as extended vector field. Indeed, ii q = (p, v) is a point 
of tangent bundle TM , then one can map it to the vector v G Tp{M). Now we can 
calculate modulus of velocity vector v and scalar product of vectors v and F. In 
terms of force field F differential equations (1.3) are written as: 



x'' = v'', 



^tv'^^F'^ix',... ,x^,v',... ,«"). 



(1.7) 



Here Vtw*^ are components of vector Vtv, where Vt is a covariant derivative with 
respect to time variable t along trajectory: 

n n 

1=1 j=l 

Definition 1.4. Newtonian dynamical system in smooth Riemannian manifold M 
is a dynamical system determined by some extended vector field F in M. 
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Note that Newtonian dynamical systems (1.7) in Riemannian manifolds are not 
purely artificial objects obtained as mathematical generalizations of the equation 
(1.1). As shown in Chapter II of thesis [1], they arise in describing constrained me- 
chanical systems with holonomic constraints. Riemannian metric in configuration 
space of such systems is given by quadratic form of kinetic energy. 

2. Extended tensor fields. 

Extended tensor fields are defined in a similar way as extended vector fields in 
definition 1.3. Let's denote by Tg{p,M) the following tensor product: 



T^{p, M) = Tp{M) eg) ... eg) Tp{M) (g t;(m) (» . . . (g r;(M) 

s times 

Tensor product T^(p, M) is the space of tensors of type (r, s) at the point p G M. 

Definition 2.1. Extended tensor field X of type {r,s) in M is a tensor- valued 
function that maps each point q ^ G C TM to a tensor of the space Tg{p,M), 
where p = Tr(q). Subset G of TM is a domain of extended tensor field X. If 
G = TM, then X is called global extended tensor field in M. 

Traditional tensor fields of type (r, s) in M are sections of tensor bundle T^M. 
Extended tensor fields of type (r, s) are sections of pull-back tensor bundle tt* (T^M) 
induced by the map of canonical projection tt : TM — > M. Below we recall some 
facts concerning extended tensor fields. Detailed explanation of the theory of such 
fields can be found in Chapters II, III, and IV of thesis [1]. 

The most important fact of the theory of extended tensor fields in Riemannian 
manifolds is the presence of two covariant differentiations 

V: r:(M)^T;+i(M), V: t;(m) ^ t;+i(m). 

First covariant differentiation V is called spatial differentiation or spatial gradient. 
In local coordinates it is represented by formula 

r n s n v / 



I \ ^ \ ^ pife Y'ii...ak...ir _ \ ^ \ ^ ybk vii ir- 

"•" Z^ Z^ <l°'h ^ji js Z^ Z^ qjk 3i---bk---js' 

k—1 afc— 1 k—1 h^ — 1 

Second covariant differentiation V is given by much more simple formula: 

v,^;-;;}: - -f^- (2.2) 

It is called velocity differentiation or velocity gradient. Velocity gradient V is de- 
fined in arbitrary smooth manifold. Unlike V, it doesn't require the presence of 
Riemannian metric in the manifold. 
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3. COVARIANT REPRESENTATION OF EXTENDED TENSOR FIELDS. 

Note, that if we replace tangent bundle TM by cotangent bundle T*M, we obtain 
another definition of extended tensor fields in M. 

Definition 3.1. Extended tensor field X of type {r,s) in M is a tensor-valued 
function that maps each point q ^ G C T*M to a tensor of the space TJ (p, M) , 
where p — 7r(g). Subset G of T*M is a domain of extended tensor field X. If 
G = T*M , then X is called global extended tensor field in M . 

In the case of arbitrary smooth manifold M definitions 2.1 and 3.1 lead to differ- 
ent theories. But for Riemannian manifold M tangent bundle TM and cotangent 
bundle T*M are bound with each other by duality maps 

g: TM -> T*M, g-^ : T*M -> TM. (3.1) 

In local coordinates duality maps (3.1) are represented as index lowering and index 
raising procedures applied to the components of velocity vector v: 



c—1 a=l 

Due to duality maps (3.1) two objects introduced by definitions 2.1 and 2.2 are 
the same in essential. We call them contravariant and covariant representations of 
extended tensor field X. If we take covariant representation of X, then formula 
(2.1) for spatial covariant differentiation V is rewritten as 

a=l b=l " 

r n s n 

I \ ^ \ ' -nik Y-il---a.k---ir __ \ ' \ ' pfcfc vil 

+ Z^ /^^qak^h Os Z^ 2^^q]k^3i...hk- 

fe=lafc = l k=lbk = l 

Formula (2.2) for velocity gradient V now is written as follows 

Si. 

dvk 



(3.2) 



v,^J;::}: = E5.^-^- (3-3) 



In order to make formulas (2.2) and (3.3) more similar to each other we raise index 
q in (3.3). As a result we get the following formula for V: 

aj^ii...ir 

V^Xl-';= ''■■■'' . (3.4) 

4. Differentiation along curves. 

Let p = p{t) be some parametric curve (e. g. the trajectory of Newtonian dy- 
namical system (1.7)). Suppose that at each point p{t) of this curve some tensor 
X = X(i) of type (r, s) is given. If X(t) is smooth function of t, then one can 
differentiate it with respect to parameter t along the curve. This is done by means 
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of covariant derivative V^. As a result we get another tensor- valued function V4X 
on the curve. Its components are given by the following well-known formula: 



Jl-.-Js I V^ V^ V^ pifc Yil...a^...ir -i 

^3l-3s - fli '^ 2^ 2^ 2^^ q^k^ji 3s-^ 

q—1 k—1 Qfc — 1 



(4.1) 



z^ A^ 2-^ lik ^ji...bk...j, 

q=l k=l fcfe = l 



X^. 



Formula (1.8) is a special case of formula (4.1) with X = v(t), and with time 
derivatives x'^ being replaced by v'^ . 

Now suppose again that some curve p ~ p{t) in M is given. Its tangent vector v 
with components x^, ... , i" is vector- valued function of parameter t. Taking pairs 
q = (p, v), where p = p{t) and v = v(i), we construct a parametric curve q — q{t) 
in TM . This curve is called natural lift of initial curve p = p(t). 

Suppose that X is some extended tensor field of type (r, s) in M. According to 
the definition 2.1, it is tensor- valued function with argument q € TM. Substituting 
q = q{t) into the argument of extended tensor field X((7), we get tensor- valued 
function X(t). li q = q{t) is natural lift of curve p — p{t), then tensor-function X(t) 
is called natural restriction of extended tensor field X to the curve p — p{t). Let's 
apply Vt to X(t). As a result for components of tensor field VtX we get 

n n 

fe=i fc=i 

One can easily write formula (4.2) in coordinate free form. Here it is: 

VtX = C(VX«)v) + C(VX®Vtv). (4.3) 

Note that V and V in right hand sides of formulas (4.2) and (4.3) are spatial and 
velocity gradients respectively, while C is the operation of contraction. 

5. Lagrangian dynamical systems. 

Lagrangian dynamical system is a special case of Newtonian dynamical system. 
The equations of dynamics (1.3) for them are given in implicit form by equations 

''"^-1^-0, fc.l,...,n. (5.1) 



dt dx^ dx'^ 

Here L ~ L{x^,... ,x"-,x^,... ,x") is Lagrange function. Differential equations 
(5.1) are known as Euler- Lagrange equations. Differentiating composite function, 
we can rewrite Euler-Lagrange equations as follows 



^ d'L ., , ^ d'L . dL 

= 1 " ■ " " s=l 



Ed L ^ oL _ 



^-^ dv^ dv'^ ^-^ dv^ dx^ dx^ 



It's clear that L is a scalar function in TM. In other words, it is extended scalar 
field. Therefore we can rewrite (5.2) in terms of covariant differentiations deter- 
mined by formulas (2.1) and (2.2). Let's use (1.8) for to express v^ in (5.2) through 
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covariant derivative Viu". Then let's use formula (2.1) applied to scalar field L for 
to express partial derivative dL/dx^ through spatial gradient 'SI kL. This yields 

n q2j^ n n n q2 j^ 

1^ Qyk QyS ■ ^f"' -l^l^l^ ^h Q^k QyS ^* '"^ + 
s — 1 i—1 7 — 1 S— 1 

(5.3) 

'^ dv^ dx" " ^L^A^ 1^3 Qy3 

s—1 S—1 J— 1 

Gathering second, third, and fifth terms in (5.3) and using formulas (2.1) and (2.2), 
we can transform (5.3) to the foUowing form: 

n n 

Y, VsVfcL • Vtv' + Y, ^s^kL ■ v' - VkL - 0. (5.4) 



Now, if we recall formula (4.2), we can further simplify our equations (5.4): 

Vt(Vfci)-Vfci = 0. (5.5) 

This form of Euler-Lagrange equations is quite similar to initial one. But now these 
equations are written in terms of covariant derivatives (2.1) and (2.2). 

6. Legendre transformation. 

In what case Lagrangian dynamical system determined by Euler-Lagrange equa- 
tions (5.5) can be written in Newtonian form (1.7)? The answer to this question 
depend on the value of determinant of the matrix A with the following components: 

dL 
^^.=V.V,^.^-^. (6.1) 

If det A ^ then, using (5.4), we can express Vtw* in explicit form, thus obtain- 
ing expression for the force field of corresponding Newtonian dynamical system. 
Lagrangian dynamical system for which the condition 

dot A 7^0 (6.2) 

is fulfilled is called regular. Now suppose that L is Lagrange function for regular 
Lagrangian dynamical system (5.5). Then let's denote p — \7L. It is clear that p 
is an extended covector field with components 

Pk=SIkL^ —^, 1, ...,n. (6.3) 

Extended covector field p with components (6.3) is used to define nonlinear map 
A: TM -^ T*M. Indeed, if g == {p,v) is a point of TM, then, taking p = p(gr), we 
can construct another pair q = {p, p) being a point of T*M: 

X{q)=q={7r{q),p{q)). (6.4) 
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Nonlinear map A : TM -^ T*M defined by formula (6.4) is called Legendre trans- 
formation. The above condition (6.2) provides local invertibility of Legendre trans- 
formation. Traditionally matrix (6.1) is assumed to be a positive matrix (see [2]): 



^>0. 



(6.5) 



This means that A is a matrix of positive quadratic form. Under the condition (6.5) 
Legendre transformation (6.4) is globally invertible, i.e. it is nonlinear bijective map 
Tp{M) -^ T*{M) at each point p £ M. The whole set of maps binding tangent 
and cotangent bundles (including linear maps (3.1)) is shown on diagram below: 



TM 



TM 




T*M 



T*M 



In section 3 above we agreed to make no difference between covariant and con- 
travariant representations of extended tensor fields. This means that we consider 
X and Xog as two forms of the same object. Differentiations V and V for these 
two representations of X are defined by formulas (2.1), (2.2), (3.2), and (3.4). Dif- 
ferentiations V and V satisfy the following identities: 



V(Xog) = (VX)og, 



V(Xog) = (VX)og. 



(6.6) 



Presence of nonlinear maps A and A^^ on diagram above increases the number of 
representations of extended tensor field X. If Y = X o A~^, then we say that Y is 
p-representation or momentum representation for X, while X is v-representation or 
velocity representation for Y. For instance, VL is a v-representation for covector 
field p, which is called the field of momentum. In p-representation components of 
covector field p are treated as independent variables pi, ... , p„. 

Legendre transformation A does not commute with differentiations V and V. 
Unlike (6.6), here we have the following equalities: 

V(YoA) =C(VYoA®VVi), (6.7) 

V(YoA) = VYoA + C(VYoA®VVi). (6.8) 

In local coordinates the equalities (6.7) and (6.8) are written as follows: 

n 

'^rX'^::X = E ^-^fc^ • ^'^.::t ' (6-9) 

n 

fc=i 
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Here X = Y » A. Similar to (4.2), formulas (6.9) and (6.10) express the rule of 
differentiation for composite functions. They are proved by direct calculations. 
Further let's define the following extended scalar field in v-representation: 

n 

h^^v^ -VkL- L. (6.11) 

fc=i 

Then let's take composition of h with Legendrc map A^^: 

H = ho\-\ (6.12) 

As a result we get another extended scalar field H. It is called Hamilton function. 
Hamilton function is traditionally used in covariant p-representation. This means 
that its argument is a point q = {p, p) of cotangent bundle T*M: 

H = H{x^,.. . ,a:", pi,... ,Pn)- 

Due to (6.12) we have h = H o X. First let's calculate Vih directly, using formula 
(6.11), then let's apply formula (6.9) for the same purpose. This yields 



fc=i fc=i fe=i 

n 
k=l 

Comparing these two formulas for V^/i and taking into account that matrix A with 
components (6.1) is non-degenerate, we obtain 

v''=W''H, k^l,...,n. (6.13) 

Formula (6.13) is analogous to (6.3), it yields an explicit expression for inverse 
Legendre map A^^ in local coordinates. Moreover, this formula means that VH is 
a p-representation for vector field of velocity v. 

Matrix A with components (6.1) is Jacobi matrix for Legendre map A. Using 
(6.13), we can calculae components of Jacobi matrix B for inverse Legendre map: 

BH = V^V^iJ. (6.14) 

Matrix B with components (6.14) inherits properties of matrix A, i.e. deti? ^ 0, 
and if A > 0, then B > 0. 

Let's denote p-representation of L by L Then let's transform (6.11) to p-re- 
presentation. Note that p-representation for v is VH and p-representation for VL 
is p. Therefore, combining (6.3), (6.13), (6.11), and (6.12), we derive 

n 

l^^Pk-^'^H-H. (6.15) 

fc=i 
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Since / by definition is p-representation of L, we can write the equality 

L^loX. (6.16) 

Formulas (6.15) and (6.16) are quite similar to (6.11) and (6.12). This reflects 
symmetry of direct and inverse Legendre maps. 

Let's calculate Vih. First let's do it directly, using formula (6.11). Then let's 
apply formula (6.10). As a result we get 



fc=i 

n 

\/^h = V^H + Y^ Vr^kL ■ \/^ H . 

fc=l 

Comparing these two formulas for V^/i, let's take into account (6.13). This yields 

\7^H = -V,L. (6.17) 

In coordinate free form the equality (6.17) is written as 

VH ^-{VL)oX-\ (6.18) 

Formula (6.18) also reflects the symmetry of direct and inverse Legendre maps. 

7. Hamiltonian dynamical systems. 

Legendre transformation is used in order to write Lagrangian dynamical system 
(5.5) in Hamiltonian form. Note that v'' is time derivative of x''. Therefore the 
above equations (6.13) are the equations of dynamics by themselves: 

x'= = V'=if(x\...,a;",pi,... ,p„). (7.1) 

But they are not complete. In order to complete these equations (7.1) we are to 
calculate time derivatives forpi, . . . , p„. Let's do it using formulas (6.3) and (4.2): 

n n 

^tPk = Y^ V.VfeL • Vtv' + Y ^s^kL ■ v'. (7.2) 

s=l s=l 

Comparing (7.2) and (5.4) and using formula (6.17), we obtain 

VtPk = -^kH{x\ . . . , x", pi, . . . ,pn). (7.3) 

Now we see that the equations (7.1) and (7.3) form complete system of ODE's. 
They are called Hamilton equations. We gather them into a system: 

x'^ = V'^il, VtPk = -VkH. (7.4) 

Hamiltonian dynamical system given by the equations (7.4) is called regular if 
matrix B with components (6.14) is non-degenerate. Each regular Hamiltonian 
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dynamical system (7.4) is locally equivalent to some regular Lagrangian dynamical 
system (5.5), and vice versa, each regular Lagrangian dynamical system (5.5) is 
locally equivalent to some Hamiltonian dynamical system (7.4). If Legendre map 
defined by (6.13) is bijective, then this equivalence is global. 

8. FiBERWISE SPHERICALLY SYMMETRIC LAGRANGIANS. 

Extended tensor field X is called fiberwise spherically symmetric if its compo- 
nents depend only on modulus of velocity vector: 

This means that X^'''^ is spherically symmetric function within each fiber of tan- 
gent bundle TM for each fixed point p € M. Such fields were considered in Chap- 
ter VII of thesis [1]. Now suppose that Lagrange function L of some Lagrangian 
dynamical system (5.5) is fiberwise spherically symmetric scalar field. Let's write 
this field in Newtonian form and let's calculate its force field F. Components of F 
should be obtained from the equations (5.4) written as follows: 

n n 

J2^s^kL-F'+J2^s^kL-v' -VkL = 0. (8.1) 

s=l s=l 

Let's calculate covariant derivatives V^V^L and VgVfcL, assuming L to be fiberwise 
spherically symmetric. For first order derivative \7kL we have 

VfcL-L'-^, (8.2) 

|v| 

Here by L' we denote partial derivative of the function L{x^, ... ,x",v) with respect 
to its last argument v, which is interpreted as modulus of velocity vector v. Let's 
apply covariant derivatives Vg and Vg to (8.2). This yields 

VgVfcL = VsL' ■ ^, (8.3) 

|v| 

VgVfeL = L ■ -—-- -f- -p-r • 5sfc - -p-p^ . (8.4) 



Note that (8.4) are components of matrix A (see (6.2) above). In order to invert 
this matrix let's consider two operators of orthogonal projection Q and P: 

First of them is a projector to the direction of velocity vector v, second is a projector 
to hyperplane perpendicular to v. Projection operators Q and P with components 
(8.5) are complementary to each other, this means that 

Q + P = l and QoP = PoQ = 0. (8.6) 
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Comparing (8.4) with (8.5), we find that 

Ask = 'VsVkL^L" -Qsk + ^^-Psk- (8.7) 

Matrix A with components (8.7) is non-degenerate if and only if L' 7^ and L" ^ 
simultaneously. In this case matrix B = A^^ has the following components: 

Now, combining (8.1), (8.3), and (8.8), we derive formula for components of F: 
This formula (8.9) is quite similar to the following one: 

s—l \ I I / 

Here W is some fiberwise spherically symmetric scalar field with W 7^ 0. Force 
fields of the form (8.10) arise in the theory of Newtonian dynamical systems admit- 
ting normal shift (see Chapter VII of thesis [1]). Let's find in which case formulas 
(8.9) and (8.10) do coincide. This occurs if the following equations hold: 

L" |v| • L" L' W ' L' W ■ ^ ' 

We consider nontrivial case, when VVK 7^ and VL 7^ 0. In this case equations 
(8.11) mean that spatial gradients of scalar fields L and L' are coUinear. This 
occurs if and only if L' = f{L, v), where / = f{u, v) is some smooth function of two 
variables. This means that we deal with a class of functions L = L{x^, . . . ,x",v), 
each of which is a solution of partial differential equation 

with some particular function / = f{u,v) 7^ 0. Note that we should not solve 
the equation (8.12) for particular function /. We should describe the whole set of 
solutions for all equations of the form (8.12). This is done by formula 

L^(3{C{x^,... ,x"),v). (8.13) 

Here (3 = (3{u,v) is a smooth function of two variables with /?(, 7^ and /3"^ 7^ 0, 
while C = C(a;^, . . . , x") is a function of spatial variables only, i. e. this is traditional 
(not extended) scalar field in M. Let's substitute (8.13) into the equation 

VsL' VsL VsL 
L" ^ Ivl • L" ^ IT' 
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derived from (8.11). This leads to the foUowing differential equation for (3{u,v): 

B" 3' 0' 

B" V B" ~ 0' ' 

This equation can be transformed so that it can be further integrated: 



0" 1 0" 



(8.14) 



Integrating (8.14) with respect to w, we obtain 

log(/3;) = log(w)+log(/3;)+log(c), where c^c{u). (8.15) 

Note that, ultimately, in formula (8.13) for L we substitute u ~ C(x^,... ,x"), 
where C(a;^, . . . , x") is arbitrary smooth function. This means that varying func- 
tion c{u) in (8.15), we do not change class of Lagrange functions L. Let's choose 
c{u) = 1/u for the sake of further convenience. Then we get 

u-0i^v0',. 

This equation is explicitly integrable. Its general solution is determined by one 
arbitrary smooth function of one variable <j> — (t>{z): 

0{u,v)^(l){u-v). (8.16) 

Substituting (8.16) into (8.13) and further into (8.9), we obtain 



Fr^^\^^-{2v^Vr-\^\^-5l). (8.17) 

s— 1 

If we substitute C — e^-^, where / = /(x^, . . . ,x"'), we can rewrite (8.17) as 

n 
Fr = -|v|2 . V,/ + ^ 2 ■ {Vsfv') ■ Vr (8.18) 



Newtonian dynamical system (1.7) with force field F given by formula (8.18) coin- 
cides with geodesic flow of metric g — e^^-^ • g, which is conformally equivalent to 
basic metric g of Riemannian manifold M. Thus we have proved a theorem. 

Theorem 8.1. Newtonian dynamical system with force field of the form (8.10) 
possess Lagrangian structure with fiber wise spherically symmetric Lagrange function 
L if and only if its force field is given by formula (8.18), which is special case of 
formula (8.10) with W = v ■ e^^ , where f = f{x^, . . . , x"). 
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9. Inverse problem of Lagrangian dynamics. 

Theorem 8.1 is not an ultimate result concerning Lagrangian structures of New- 
tonian dynamical systems admitting normal shift. First reason is that formula 
(8.10) does not cover general case (see Chapter VII of thesis [1]). General for- 
mula for the force field of Newtonian dynamical system admitting normal shift of 
hypersurfaces in Riemannian manifold M with dim M ^ 3 looks like 

where h = h{w) is an arbitrary smooth function of one variable. And second reason 
is that in theorem 8.1 we restrict ourselves to the case of fiberwise spherically 
symmetric Lagrange functions. 

In order to study general case we should substitute (9.1) into (8.1) and consider 
(8.1) as a system of PDE's for unknown Lagrange function. Problem of determin- 
ing whether the Newtonian dynamical system with a given force field F admits 
Lagrangian structure (and finding Lagrange function if it admits) is known as in- 
verse problem of Lagrangian dynamics. As known to me, this problem is not solved 
in general case (see more details and references in [3-7]). Even for special force 
fields given by explicit formula (9.1) it remains unsolved. Solving this problem for 
force field (9.1) is very important since it would open a way for applying theory 
from [1] to the description of wave front dynamics and to some related problems 
arising in analysis of partial differential equations (see [8] and [9] ) . 

10. Resume. 

Concept of extended tensor field arisen in [10] and [11], and used in [1] for 
describing Newtonian dynamical systems is applicable to Lagrangian and Hamil- 
tonian dynamical systems in Riemannian manifolds as well. It gives a method (or 
a language) for describing these systems in terms of their configuration space M 
instead of using geometric structures in tangent bundle TM (exception is Hamilton- 
ian dynamical systems in abstract simplectic manifolds, when one cannot separate 
configuration space within phase space of dynamical system). As an example of 
applying suggested method I consider inverse problem of Lagrangian dynamics for 
Newtonian dynamical systems admitting normal shift, and I give partial solution 
of this problem in class of fiberwise spherically symmetric Lagrange functions. 
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